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ABSTRACT 

The nearby star a Oph (Ras Alhague) is a rapidly rotating A5IV star spinning at ^89% of its 
breakup velocity. This system has been imaged extensively by interferometric techniques, giving a 
precise geometric model of the star's oblateness and the resulting temperature variation on the stellar 
surface. Fortuitously, a Oph has a previously known stellar companion, and characterization of the 
orbit provides an independent, dynamically-based check of both the host star and the companion 
mass. Such measurements are crucial to constrain models of such rapidly rotating stars. In this 
study, we combine eight years of Adaptive Optics imaging data from the Palomar, AEOS, and CFHT 
telescopes to derive an improved, astrometric characterization of the companion orbit. We also use 
photometry from these observations to derive a model-based estimate of the companion mass. A fit 
was performed on the photocenter motion of this system to extract a component mass ratio. We find 
masses of 2.40+037 M Q and O.85+004 M Q for a Oph A and a Oph B, respectively. Previous orbital 
studies of this system found a mass too high for this system, inconsistent with stellar evolutionary 
calculations. Our measurements of the host star mass are more consistent with these evolutionary 
calculations, but with slightly higher uncertainties. In addition to the dynamically-derived masses, 
we use IJHK photometry to derive a model-based mass for a Oph B, of 0.77±0.05 M Q marginally 
consistent with the dynamical masses derived from our orbit. Our model fits predict a periastron 
passage on 2012 April 19, with the two components having a 50 mas separation from March to May 
2012. A modest amount of interferometric and radial velocity data during this period could provide 
a mass determination of this star at the few percent level. 

Subject headings: instrumentation: adaptive optics — stars: individual (HIP86032) 



1. INTRODUCTION 



Although the binar y prop er ties of solar-type stars 
(iDuquennov fc Mavorl 119911: iMetchev fc Hillenbrandl 
2009 ) and lower mass late -type stars ()Fischer fe Marcvi 
1992; IReid fe Gizisl Il997h . are becoming increasingly 



I Department of Astronomy, California Institute of Technol- 
ogy, 1200 E. California Blvd, MC 249-17, Pasadena, CA 91125 

University of Michigan, Astronomy Department, 941 Denni- 
son Bldg, Ann Arbor, MI 48109-1090 

3 Astrophysics Department, American Museum of Natural 
History, Central Park West at 79th Street, New York, NY 10024 

4 Jet Propulsion Laboratory, California Institute of Technol- 
ogy, 4800 Oak Grove Dr., Pasadena CA 91109 

5 Department of Astronomy, Columbia University, 550 West 
120th Street, New York, NY 10027 

6 Institute of Astronomy, University of Cambridge, Madingley 
Road, Cambridge CB3 OH A, UK 

7 School of Physics, University of Sydney, NSW 2006, Aus- 
tralia 

8 National Astronomical Observatory of Japan, Subaru Tele- 
scope, Hilo, HI, 96720. 

9 CFHT Corp., 65-1238 Mamalahoa Hwy., Kamuela, HI 96743 

10 STScI, 3700 San Martin Drive, Baltimore, MD 21218 

II Stony Brook University 

12 NSF Center for Adaptive Optics. 

13 Caltech Optical Observatories, California Institute of Tech- 
nology, Pasadena, CA 91125 

14 NASA Exoplanet Science Institute, California Institute of 
Technology, Pasadena, CA 91125 

15 Department of Astronomy, Cornell Univ., Ithaca, NY 14853 

16 UCLA Department of Astronomy 

17 Sagan Fellow 

18 NSF Postdoctoral Fellow 



well understood, studies of the higher end of the mass 
spectrum are crucial for a broader understanding of stel- 
lar multiplicity. Indeed, recent simulations suggest that 
A stars may possess a greater ab undance of both stel lar 
and planetary mass companions (Kratt er et al.ll2010L J. 
Crepp, private communication). However, the multiplic- 
ity characteristics of these more massive stars has only 
recen t ly started to be su r veyed (|Shatskv fe Tokovininl 
I2002t IKouwenhoven et al.l I2005D . and more extensive 
surveys are underway to place limits on the nature of 
massive-star multiplicit y. In recent years, high c o ntrast 
imaging techniques ([Oppenheimer fe Hinklevi 120091 : 
lAbsil fc MawetJl2009D such as adaptive optics (hereafter 
"AO" ) have matured significantly, allowing these studies 
to proceed. The nearby (14.68 pc) A5IV st ar a Oph 
(Ras Alhague) has a well known companion (|WagmarJ 
ll946tlLippincott fc Wagmanll 19661: lGatewoodll2005[ ) with 
a 8.62-year period, well established o yer several decade s 
of monitoring and first resolved by iMcCarthvl (|1983ft . 
But a fuller characterization of the companion has not 
been possible, and, aside from the peri od, the orbital 
parameters have been loosely constrained dKamper et al.1 
[l^lAugensen fc Heintdll992t [Gatewoodll2005l) . ~ 

Characterization of the binary nature of A-star sys- 
tems becomes doubly interesting when the host stars are 
also rapidly rotating. Recent interferometric imaging of 
several rapidly rotating A-stars have revealed imagin g of 
their surfaces (jMonnier et al.ll2007t IZhao et al.ll2009f ). a 
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TABLE 1 
Table of Astrometry 



epoch 


MJD 


r (mas) 


(7 r 


PA 


CPA 


x, y (mas) 


Band 


Telescope (Instrument) 


Rcf a 


1 


2445157.9 










-380Y±50 


K,L 


3.8m KPNO Mayall 


1 


2 


2445244.8 










-650X±50 


K,L 


3.8m KPNO Mayall 


1 


3 


2445424.2 










-340Y±50 


K,L 


3.8m KPNO Mayal 


1 


A 
'1 


O A A KQCiA fl 

2440oU4.U 










-430Y±oO 


K , L 


o.om ixr^lNU Mayall 


1 


5 


2447045.8 










-400X±50 


K, L 


3.8m KPNO Mayall 


1 


6 


2451304.5 


770 


± 40 


243.7 


± 3.0° 




2.15/im 


3.6m La Silla (SHARP II) 


2 


7 


2452473.8 


470 


± 20 


233.2 


± 2.4° 




I 


AEOS (Vislm) 


3 


8 


2452760.0 


<275 


± 3 


n/a 






I 


AEOS (Vislm) 


3,4 


9 


2453168.5 


303 


± 33 


253.0 


± 6.3° 




H 


AEOS (Lyot Project) 


5 


10 


2454253.9 


776.5 


± 2.1 


244.6 


± 0.4° 




J,H,K 


Palomar (PHARO) 


6 


11 


2454263.5 


765 


± 20 


243.7 


± 1.5° 




H 


AEOS (Lyot Project) 


5 


12 


2454637.5 


787.8 


± 2.8 


240.6 


± 0.4° 




K, Br 7 


Palomar (PHARO) 


6 


13 


2454657.5 


790 


± 20 


239.5 


± 1.4° 




J,H 


Palomar (Project 1640) 


7 


11 


2454963.8 


756.6 


± 7.5 


239.3 


± 1.2° 




Pap (1.28/itn) 


CFHT (PUEO) 


8 


15 


2455002.5 


748.5 


± 2.7 


238.3 


± 0.4° 




K 


Palomar (PHARO) 


6 



1 References: 



: 1 ) IMcCarthvl i 19831) and IKamper et ahl 1119891) provide one-dim ensional speckle measurements 2) 
Boccaletti eTaLl 120011 ), 3) IRoberts fc Nev manl ||2002 | 1, 41 Upper lim it to orbit, 5) Sivaramakrishnan et al. (2007]), 6) 
Havward et al.l 120011 )7 7) IHinklev et al.l 12008) . 8) IRigaut et all (WM) . 



Oph, a prototypical rapid rotator, is rotating at ~ayyo 
of its breakup velocity. Interferometric surface imaging 
of this star clearly shows its pronounced oblateness: the 
star's equatorial radius is 20% larger than its polar ra- 
dius, con^s_rjondhi£_to a ^1840K temperature differen- 
tial (|Zhao et al.ll2009D . This 2-3 M star has also been 
the target of extensive asteroseismic monitoring carried 
out by the MOST satellite, discovering ~50 pulsational 
modes (J. D. Monnier, private communication). 

Derivation of the component masses for this system is 
a critical check against investigations of the host star's 
rapid rotation and asteroseismology. Discrepancies be- 
tween the dynamically-derived mass of the host star and 
that allowed by models of rapid rotators, may illuminate 
potential complications with rotator models, e.g. those 
that do not take into account differential rotation of the 
star. For all these reasons, deriving a more accurate or- 
bit is crucial for future studies to reconcile the orbital 
dynamics of the companion star with the physics dic- 
tating the rotation of the host star. In this paper, we 
assemble data from a number of high-contrast imaging 
programs over the past eight years. We take advantage 
of these recent measurements along with archival data to 
further constrain the orbit of the companion, as well as 
the system mass. 

2. ASTROMETRIC OBSERVATIONS 

Table Q] presents astrometric measurements of a Oph 
obtained using several instruments. Epochs one through 
five are p reviously measured o ne-dimensional speckle 
data from IKamper et al.l (JT989). For the sixth epoch 
of observations, we adopt the radial separation (770± 
40 mas) and position angle (2 4 3. 7± 3.0 °) measurements 
report ed in iBoccaletti et al.l (|2001l ). iBoccaletti et al.l 
(2001) mistakenly reported t he position angl e in de- 
grees measured West of North. iGatewoodl (|2005t ) noticed 
the problem, but mistakenly assumed it was a quadrant 
problem and changed the value by 180°. We use the cor- 
rect measurement position angle of 243° (A. Boccaletti, 
private communication). 

Epochs seven and eight are comprised of visible AO 
measurements obtained from the 3.63m AEOS tele- 
scope on Haleakala, Hawaii {R oberts fc Nevmanl 




Astrometric Orbit for 
Alpha Oph 

200 -200 -400 -600 -800 

Milliarcseconds 

Fig. 1. — The visual orbit of the a Oph system. The solid line is 
the best fit solution while the grey band shows the allowed orbits 
based on Monte Carlo sampling of the parameter uncertainties. 
The primary is located at (0,0) and the astrometric data from 
Table [l] arc plotted with the positional uncertainties. East is left, 
and North is up. 



2002; IRoberts et all l2005h . Epochs nine and 
eleven were obtained using " The Lyot P r oject" 
dSivaramakrishnan et "all 120071: IHinklev et "all 120071: 



iLeconte et al.l 20101). a d i ffraction-limited classical Lyot 
coronagraph (|Lvotlll939t ISivaramakrishnan et al.l 120011) 
working in the infrared and recently decommissioned at 
AEOS. 

The thirteenth epoch of observations of the a Oph sys- 
tem was obtained using a recently commissioned coro- 
nagraph integrated with an integral field spectrograph 
(IFS) spanning the J and 77-bands (1.10/im - 1.75/im) 
on the 200-in Hale Telescope at Palomar Observatory 
(Hinkley et al. 2008). This instrument package, called 
"Project 1640" is mounted on the Pa lomar AO system 
(]Dekanv et alj fl998t iTrov et al.l l2000ft . and is a dedi- 
cated high contrast imaging instrument providing low 
resolution spectra (A/AA)~30 — 60 with a len slet-based 
IFS and Apodized-Pupil Lyot coronagraph (jSoummerl 
2005), an improvement of the clas sical Lyot corona- 
graph (jSivaram akrishnan et al. 2001). Epoch 14 obser- 
vations wei^_jnadeat_the Canada-France-Hawaii Tele- 
scope (|Rigaut et al.l 119981) . The observations compris- 
ing epochs 10, 12, and 15 were obtained with the 
Palomar AO system and the PHARO infrared camera 
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TABLE 2 

Binary Orbit and Derived Properties 



Parameter 


Value 


From Visual Orbit Only: 




Semi-major axis a (mas) 


427Ij3 


Eccentricity 


0.92 ± 0.03 


Inclination (degs) 


125±« 


uj (degs) 


162 ±14 


O fHeixs'l 
1l L 1 LLC 1 


232±9 


Period (days) 


3148.4 (fixed) 


Epoch of pcriastron (JD) 


2452888 ± 53 


From Astrometric Photocenter 




Motion (Using Visual Orbit Above): 




Mass Ratio (Primary/Secondary): 


2.76 tf 7 


Parallax (mas) 


69.1 (fixed) 


Fixed proper motion! E,N (mas yr ) 


1 OQ Q 997 n 


System Mass (M ): 


3-251 H 


a Oph A (Mq): 


2.40+3 7 


a Oph B (Mq): 


0.85±^ 


TABLE 3 
Apparent Magnitudes for the a 


Oph system 



Member 


J 


J 


H 


K 


a Oph A a 


1.90 


1.83 ±0.28 


1.72 ±0.18 


1.68 ±0.21 


a Oph B_ 


6.72 ± 0.10 


6.02 ± 0.310 


5.44 ±0.206 


5.25 ±0.236 



a 7-band photometry: IMonet et all p003ft , JHK: ICutri et all 1 120031 
b /-band photometry was obtained using the AEOS Vislm instrument 
and JHK photometry with the PHARO infrared camera at Palomar. 

(|Havward et alJfeOOlf) . Astrometric points from epochs 
six through fifteen are shown in Figure [1] and represen- 
tative examples of data from epochs 9, 11, and 13 are 
shown in Figure [2] 

3. ORBITAL CHARACTERIZATION 

We fitted an orbit to the ast rometric data in Tab l e IT ] 
including the speckle data from lKamper et all (|198 9) Q 
We used both a Levenberg-Marquardt method (using 
IDL procedure MPFIT developed by C MarkwardtET 
and also the IDL routine AMOEBA (jPress et al.lll992 
arriving at the same global solution. Comparing our pre- 
dicted astrometry with the data, we find a reduced % 2 
of 0.7. The best-fit orbital elements can be found in Ta- 
ble [H where we have fixed the quad r ant o f lo using the 
radial velocity data of iKamper et al.l ()1989f ). In order to 
calculate errors for the orbital elements, we carried out a 
Monte Carlo simulation, allowing all parameters to vary 
in the fitting except for the period, which we fixed to 
8.62 years. The results of the error analysis can be found 
in Tabled 

Using the full covariant set of parameters from the 
monte carlo study, we were then able to determine the 
component mass ratio by fitting the p hotocenter motio n 
measured with the MAP instrument (jGatewoodl 120051 ) . 
In order to carry out a fit to the absolute position of 
the photocenter, we needed to constrain a few impor- 
tant parameters using independent sources. Firstly, we 

19 We needed to flip the signs on x,y from this work to be con- 
sistent with other astrometry measurements. 

20 http:/ /cow. physics. wisc.edu/~craigm/idl/idl. html 



adopted a proper motion of 123.3 mas yr 1 (East) and 
-227.0 mas yr -1 (North) from the FK5 catalog. Note 
that other proper motion catalogs with shorter time base- 
lines contain significantly different values for a Oph, pre- 
sumably having been corrupted by orbital motion of the 
primary. We also a dopte d a parallax of 69.1 mas from 
iMartin k Mignardl (|1998[ ). Note that the errors in these 
three quantities have been neglected in the following 
analysis (errors from our visual orbit dominate uncertain- 
ties in the mass ratio). Figure [3] shows these relative pho- 
tocenter offsets after subtracting out the system proper 
motion and an estimate of the center of mass. Lastly, 
we had to assume a flux ratio at the effect ive wavelength 
of the MAP experiment Gatewood (1987), similar to an 
i?-band filter. Based on multi-wavelength detections of 
the companion at /, J, H, and if-bands, we have esti- 
mated an i?-band flux ratio of 5.3 magnitudes (factor of 
130) assuming the companion is a mid/late K star (See 
Section l4~Tj) . Using the results from the two analyses, we 
constrain the masses to be 2.401^37 M and 0.85^^ 
Mq for a Oph A and a Oph B, respectively. Table [2] 
contains the results of this orbital fitting, along with the 
full set of the fit parameters. 

4. ANALYSIS 
4.1. Companion Spectral Energy Distribution 

Table [3] presents the broadband apparent magnitudes 
for a Oph B obtained with the PHARO infrared camera 
at Palomar, as well as the AEOS Vislm camera. The 
spectral energy distribution for a Oph B is shown in 
Figure [H accou nting for the sys tem's 0.834±0.024 dis- 
tance modulus ([Gatewo od 2005). The broadband mag- 
nitudes listed in Table [3] are the result of a simultane- 
ous photometric fit to each binary component which re- 
turns position angle, relative separation as well as rel- 
ative brightnesses. Figure @] shows /, J, H, Br 7 , and 
-ftT-band photometry. Also shown are broadband J and 
.ff-band photometric values obtained with the Project 
1640 IFS at Palomar Observatory. The Project 1640 ob- 
servations were calibrated using referenc e star spectra 
obtained from the IRTF s pectral library (Gushi ng et al.1 
120051 : iRavner et al.l 120091) in th e same manner as de- 
scribed in iHinklev et al.l (|2010D . However, these ob- 
servations were gathered during the first light observa- 
tions with this new instrument, and hence the calibra- 
tion dataset was obtained before a full understanding of 
an effective calibration process was mature. 

Also shown in Figure B] are theoretical broa dband I, J, 
H, an d if-band fluxes based on the models of lSiess et al.l 
(2000). The curves in Figure [4] are interpolations be- 
tween these values. At a spectral resolution of 30-60, 
late-G to early-M star spectra are very similar, prevent- 
ing robust discrimination between them. Rather than 
plotting template spectra over the data for comparison, 
we have instead chosen to use these model-based inter- 
polated mass-luminosity curves to allow discrimination 
of the best fit mass value for a Oph B. The best fit curve 
to the measured fluxes corresponds to a mass of 0.77 
Mq. Taking into account our photometric errors as well 
as the 0. 024 magnitude d istance modulus uncertainty 
quoted bv lGatewoo d (2005), we deduce an uncertainty of 
±0.05 M Q in our best fit value. This value of 0.77±0.05 
M is very similar to the 0.78±0.058 M Q reported by 
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Fig. 2. — Examples of data taken at the three epochs listed at the upper left hand corner of each panel. The left and middle images were 
taken with a near-IR camera with a coronagraph at the AEOS telescope. The right hand pane l is a single slice extracted from a data cube 
from the high contrast integral field spectrograph at Palomar Observatory (Hinkley ct al. 2008). The companion, a Oph B, is marked with 
a white arrow in all three epochs, and is shown near apastron in the third panel. 

This is analogous in practice to measuring the rela- 
tive orbital inclinations of trip le or quadruple system 
(e.g., Mutcrspaugh et ail I2008D which can reveal hints 
to th e formation and/or evolut ion of system dynamics 
(e.g., ISterzik fe Tokoviriml [2002V In order to carry out 
this comparison, we convert the spin parameters from 
iZhao et all (|2009t) into an equivalent Q. = 216.12 ± 1.23°, 
i = 92.30 ± 0.43° that would characterize a circular orbit 
in the star's equatorial plane. Note these values are af- 
fected by orbital degeneracies since we have not incorpo- 
rated any surface or orbital velocities into the estimate. 
However, because the primary spin geometry is nearly 
edge-on, these degeneracies do not matter much. We find 
that the mutual angle of inclination is $ = 36 ± 6° (or 
$ = 144 ± 6° for the opposite spin polarity). Inter- 
estingly, this lies close to the predicted peak inclination 
(40°) that arises from Koz ai cycles in the context of a n 
inner short-period binary dFabrvckv &: Tremainel f2007'). 
echoing findings of iMuterspaugh et al.l ( 20081 ). However, 
the primary a Oph A has no known inner companion 
and so this mechanism does not seem to be applicable. 
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Fig. 3. — T he relative photo center positio ns of a Oph from MAP 
observations (Gatcwood 1987) reported by Gatcwood (2005) after 
subtracting out the system proper motion and an estimate of the 
center of mass. With the orbit fit shown in Figure[l] the component 
mass ratio is obtained by fitting the photocenter motion shown 
here. 



iGatewoodl (|200l . and within 1.3<r of our 0.85± ^ M 
dynamical mass discussed in Section 03 Moreover, in 
this m ass regime (~0.7— O.8M ), Hillenbrand fcW hitc 
(|2004h have demonstrated t hat mass-luminosit y models, 
including those produced bv lSiess et al.l (|2000D . are only 
inaccurate at the few percent level. Due to the afore- 
mentioned uncertainties in the calibration of the Palo- 
mar Project 1640 photometry, only the Palomar PHARO 
data and the AEOS /-band data were used in the fit. Fig- 
ure 2] also shows a curve corresponding to the dynamical 
best-fit mass value of 0.85 M , as well as a curve corre- 
sponding to 0.65 M for reference. Comparing our de- 
rived near-infra red (J—H) and ( H— K ) colors with those 
documented in iBessell fc Br ett (1988) puts the spectral 
type of the companion between a K5V and a K7V. 

4.2. Spin-Orbit Mutual Inclinations 

Our measurement of the orbital plane of the Al- 
pha Oph binary system can be compared to the an- 
gular momentum of the rapidly-rotating primary star. 



5. DISCUSSION AND CONCLUSIONS 
We have obtained an orbital fit for the a Oph system 
and derived component masses of 2.40_; 37 and 0.85_; 04 
Mq for a Oph A, and a Oph B, respectively. Our 
estimation of the mass of a Oph A is lower th an the 
2.84±0.19 M value quoted bv IGatewoodl (|2005l ). How- 
ever, our Monte Carlo approach yields more conservative 
estimates of the uncertainties, overlapping with other 
wo rks. Specifically, th e 2.1 M Q value for a Oph A quoted 
by I Zhao et"aT1 (|2009f ) and measured using an HR dia- 
gram overlaid with isochrones, is allowed at this work's 
la level, while it is only allowed at the 3.5<r level in the 
lGatewoo"dl (|2005l ) study. Our work demonstrates that, es- 
pecially with a high eccentricity orbit, careful estimates 
of stellar parameters are essential. 

We anticipate that this dynamical estimation of the 
mass of a Oph A will aid forthcoming studies focused 
not only on modelling the rapid rotation of a Oph A, 
but also of the host star's asteroseismological properties 
(J. D. Monnier, private communication). The mass value 
derived here may help to distinguish between models that 
assume solid-body rotation of the star, or have some de- 
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2;ree of differential rotation beneath the photosphere. 



with a mid/late-K dwarf classification for the compan- 
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Fig. 4. — Broadband photometry for a Oph B from three of the 
observing programs described in this paper. Also shown are three 
theoretical broadband I, J, H, and if-band flu x curves based on 
the mass-luminosity models of Sicss ct al. (2000). Each curve is an 
interpolation between these values. The best fitting model is shown 
(solid line-0.77 Mq), as is the dynamically determined mass of 0.85 
Mq, and a 0.65 Mq curve for reference. The horizontal error bars 
reflect the width of each bandpass. 



Moreover, we have compiled the most complete set of 
photometric data yet for a Oph B from several observing 
programs to obtain an estimate of the companion mass 
through the use of theoretical models. The 0.77±0.05 
M Q estimation of the companion mass derived photo- 
metrically is marginally consistent with the 0.85^ 04 M 
value for a Oph B. Our infrared colors and the corre- 
sponding best fit mass-luminosity curve are consistent 



Finally, based on our fit, a Oph B will be passing 
through pcriastron on 2012 April 19, with an uncertainty 
on this date of ± 53 days. However, the separation be- 
tween the two objects should be ~50 mas from 2012 
March 24 to 2012 May 21. Such a separation will be 
ideal for existing state-of-the-art interferometers such as 
the CHARA Array or the Very Large Telescope Interfer- 
ometer. Intensive astrometric monitoring of this system 
and new radial velocity observations, especially during 
periastron, can help to constrain the mass of a Oph A 
to within a few percent. 

This work was performed in part under contract with 
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referee for several helpful comments. 
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